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Abstract. 

A general method is proposed for calculating a fully k-dependent, continuous, 
and causal spectral function A(k, E) within the recently introduced nonlocal 
version of the coherent-potential approximation (NLCPA). The method involves 
the combination of both periodic and anti-periodic solutions to the associated 
cluster problem and also leads to correct bulk quantities for small cluster sizes. We 
illustrate the method by investigating the Fermi surface of a two-dimensional alloy. 
Dramatically, we find a smeared electronic topological transition not predicted by 
the conventional CPA. 
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1. Introduction 

The coherent-potential approximation (CPA) yQ is a successful and widely used theory 
of dealing with disordered solids. However, the CPA does have a serious limitation 
in that it is only a single-site mean-field theory of disorder and hence unable to 
treat multi-site correlations in the ensemble of disorder configurations. However, 
this long-standing problem has recently been addressed by the development of the 
nonlocal coherent-potential approximation (NLCPA) [3], a cluster mean-field theory 
which provides systematic corrections to the CPA. It was introduced as the static 
version of the dynamical cluster approximation (DCA) |3J 0] used to describe short- 
range correlations in strongly-correlated electron systems, and has been implemented 
both for model Hamiltonians [3 El El and within first-principles density functional 
theory El E]. Notably, in contrast to earlier cluster theories JH] and supercell 
calculations, the NLCPA preserves the full translational symmetry of the underlying 
lattice. Significantly, this highly desirable attribute means the NLCPA could be used 
to study nonlocal effects on spectral properties of disordered systems. Indeed, the 
spectral function A b (k, E) for metallic alloys is of general interest because its peaks 
define the Fermi surface. For an ordered system these peaks are delta functions, 
however for a substitutionally disordered alloy they have finite heights and widths 
where £ is the coherence length at the Fermi energy Ep. The shapes of such Fermi 
surfaces can be determined from two-dimensional Angular Correlation of (Positron) 
Annihilation Radiation (2D-ACAR) experiments 11 , and, since they drive various 
ordering processes and are subject to observable topological changes p3 H3| . 
their theoretical description is an important area in alloy physics. One of the main 
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objectives of this paper is to introduce a method for calculating the spectral function 
within the NLCPA and to demonstrate its use by explicit calculations for the Fermi 
surface of a tight-binding model alloy. 

Recall that for such a tight-binding model alloy, the configurational average 
over all possible disorder configurations is represented by a translationally-invariant 
effective medium characterized by the exact self energy S(k). By mapping to a self- 
consistent single-site impurity problem, the conventional CPA neglects correlations 
beyond a single site and hence yields a single-site self-energy E that has no dependence 
on momentum k. The idea of the NLCPA is to systematically restore this neglected 
k-dependence by replacing the conventional single-site impurity with an ensemble 
of impurity cluster configurations, thus providing systematic corrections to the 
CPA as the cluster size increases. Unlike previous cluster theories, the NLCPA 
preserves translational-invariance by imposing Born- von Karman (periodic) boundary 
conditions on the cluster. However, a consequence of this is that the resulting effective 
medium is characterized by a self-energy S(K) which is a step function in k-space, 
with the number of steps centred at the points K equal to the number of sites in the 
cluster. This naturally leads to a spectral function with discontinuities at the jumps, 
which in turn leads for example to an unphysical Fermi surface. 

Furthermore, it has recently been explicitly shown [H] that the NLCPA results 
are different if anti-periodic |Hj Born- von Karman boundary conditions are imposed 
on the same cluster instead of the conventional periodic boundary conditions. The 
periodic and anti-periodic results only become equivalent as the cluster size iV c is 
increased to some critical value which depends on the model in question. In fact, 
below this critical value of iV c it has been explicitly shown that the periodic and 
anti-periodic solutions on their own are unphysical due to the discontinuities, i.e. not 
sufficient to yield a meaningful density of states |5] . 

In this paper we introduce a new method for implementing the NLCPA which 
resolves all the issues above. It involves the combination of both periodic and anti- 
periodic solutions to the NLCPA and should be used for cluster sizes smaller than the 
critical value. It yields a fully k-dependent, continuous, and causal spectral function 
j4(k, E), together with a meaningful density of states. 

The outline of this paper is as follows. First, we briefly summarise the idea of the 
NLCPA in section I2"T1 Then, in section f2~Sl the concept of periodic and anti-periodic 
Born-von Karman boundary conditions is described in detail, and the key idea of the 
new method for implementing the NLCPA is introduced. The strict criteria that the 
new method must satisfy are outlined in 12.31 and the formalism is detailed in section 
IO To test the validity of the new method, section 13.11 describes an application to 
a simple ID tight-binding model in order to compare with known results. Following 
this, we apply the method to a 2D tight- binding model in section Hi 21 and investigate 
the Fermi surface. Finally we conclude in section 

2. Formalism 

2. 1 . Brief summary of the NL CPA 

Consider a general tight-binding model binary alloy defined by a hopping amplitude 
Wij and site energies Va and Vb- In terms of the site occupation numbers which 
take on values of and 1 depending on whether the site is occupied by an A-atom or 
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a B-atom, the usual Green's function Gij(E) satisfies the equation 

J2 l(E - Vi)Sij - Wo] Ga = % (1) 

ij 

where Vi = & Va + (1 — &) and £ is the energy. We are interested in the average of 
Gij(E) over all possible disorder configurations Here it is useful to introduce the 
concept of an effective medium, and rewrite the average Green's function in terms of 
the exact self-energy Sjj which is translationally- invariant. In essence it replaces the 
substitutional disorder at each lattice site. By applying a Fourier transform to such 
translationally-invariant quantities we obtain 

G(k) = ( J B-H/(k)-S(k))- 1 (2) 

where a is the lattice constant and VF(k), S(k) are the hopping and exact self-energy 
respectively in reciprocal space. In the interests of clarity, the energy dependence of 
all Green's functions on energy will be implicitly defined as in J2J) for the remainder 
of this paper. 

In the single-site CPA the self-energy is k-independent and the purpose of the 
NLCPA is to systematically restore this neglected k-dependence. As explained in 
[21 E] the idea is to coarse-grain the function S(k) over a set of N c tiles which divide 
up the Brillouin zone of the underlying lattice. If these tiles are centred at a set 
of N c points {K„} then this results in a function E(K n ) which is a step function 
in k-space. The NLCPA is able to determine the function S(K„) by mapping to a 
self-consistent embedded cluster problem, where the cluster of N c sites has Born- 
von Karman (periodic) boundary conditions imposed. This means the real-space 
cluster sites {R/} and corresponding cluster momenta {K n } are related by the Fourier 
transform 

^££ e lK " (R '- Rj) = i 0) 

C K„ IJ 

where n = 1, ,.,N C . This relation is key to constructing an algorithm to determine 
the medium (see [2] for full details). After convergence of the algorithm, the NLCPA 
approximation to the exact average Green's function (J2J is given by 

G(k) = (£ - W(k) - E^))- 1 (4) 

where S(K n ) takes its appropriate value within each reciprocal-space tile n. Clearly 
the effect of coarse-graining the self-energy has been to reduce its range in real space 
(although physics on a longer length scale is treated at the mean-field level), however 
translational invariance is always preserved. 

Unfortunately, as mentioned in the introduction, the fact that E(K„) takes its 
appropriate value within each reciprocal-space tile in means that discontinuities 
have been introduced at the tile boundaries, leading for example to an unphysical 
Fermi surface. 

2. 2. Idea of the new method 

In order to remove the discontinuities in G(k) described above, an idea would be 
simply to interpolate the self-energy E(K„) to a smooth function S(k) in such a way 
as to preserve the number of states. However, the problem with such an approach is 
that new information (e.g. structure in the DOS) would be inferred. A second problem 
is that of preserving causality. For example previous attempts have been made (within 
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the DCA) |15j to interpolate the self-energy via schemes such as the maximum entropy 
method (MEM). However, although the MEM is guaranteed to be causal, it can only be 
used to add a finite number of points to a data distribution, after which it adds spurious 
structure. One must use a secondary interpolation to fill in the gaps left by the MEM 
e.g. Akima spline as used in the DCA ^Sj- However such spline interpolations are 
well known to overshoot the interpolation in some circumstances , and particulary 
when the data set used for the interpolation is rapidly varying the spline may violate 
causality. To overcome this, more data points are required and hence a bigger cluster 
size, but one is still not ensured a causal interpolation for all cases. Furthermore, in 
the case of first-principles alloy calculations |S1 03 we are particularly interested in 
large deviations from the CPA in k-space where the self-energy will vary rapidly since 
only small clusters are computationally viable at the present time. 

In this paper a different approach to the problem is introduced. First recall 
that in previous NLCPA calculations to date, periodic Born-von Karman boundary 
conditions have been applied to the cluster. This essentially means reducing the size 
of a conventional lattice with Born-von Karman boundary conditions to contain only 
a cluster of N c sites, so that the edges of the cluster map round to the other end in 
each dimension e.g. 

G [R x + L x , R y + L y , R z + L z ) — G (R x , R y , R z ) (5) 

for the cluster Green's function, where Li is the length of the cluster along direction 
i. Since the lattice constant is unchanged, the boundaries of the BZ will remain the 
same, however it will now contain only N c evenly spaced k points referred to as the 
set of cluster momenta {K„}. Thus the conventional lattice Fourier transform used 
in the N c — > oo limit reduces to the cluster Fourier transform given by However, 
as mentioned in the introduction, for a given cluster there are another set of cluster 
momenta {K„} satisfying J3J corresponding to anti-periodic 01 E] Born-von Karman 
boundary conditions so that quantities such as the cluster Green's function at the 
edges of the cluster map back to minus the value at the other end i.e. 

G (R x + L x , R y + L y , R z + L z ) = — G (R x ,R y , R z ) (6) 

As an example, the periodic and anti-periodic set of cluster momenta for a 2D square 
lattice with N c = 4 are shown in figure^ labelled as {K^} and {K^ p } respectively. 
Permutations of periodic and anti-periodic conditions along each axis are also possible, 
for example periodic along the x-axis and anti-periodic along the y-axis and vice- versa. 
These are shown in figure 2] for the same example and are labelled as {K^" 4 ^"} and 
{K^" P *} respectively. 

Next, recall that it has recently been explicitly shown |Bj that for a given cluster 
the NLCPA results change with respect to the choice of {K„}. Any theory with a 
choice of boundary condition must be invariant in its predictions of bulk properties [T7] 
and hence at present the NLCPA can only be perceived as a unique theory above some 
critical value of N c , which depends on the model in question. Moreover, below this 
critical value of N c it has been explicitly shown that the periodic and anti-periodic 
DOS results on their own are not meaningful due to the discontinuities jjj. The main 
idea of the method introduced in this paper is to appropriately combine the results 
obtained with the various boundary conditions to give a unique and meaningful result 
for any given cluster. 

First note that, significantly, the set of K„ points associated with each boundary 
condition is shifted compared to the others so that no overlap occurs. This can be 
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Figure 1. (colour online) Shown (left) are the sets of cluster momenta 
{K^} (denoted by asterisk) and {K^ p } (denoted by circles) centred at their 

P AP 

corresponding tiles. Shown (right) are the sets {K n x " } (denoted by plus) 
P AP 

and {K n y * } (denoted by trangles) centred at their corresponding tiles. Note 
that r is situated at (0, 0), X is situated at (0, ± — ) and (± — ,0), and M is situated 
at(±£,±£). 



clearly seen in figure^for the 2D square lattice example. Next, let us label the Green's 
function (given by Q) obtained using a particular boundary condition by Gj(k) and 
its corresponding set of cluster momenta by {K n }i, where i = P, AP, P x AP y P z .. 
etc. Note that on account of the coarse-graining procedure outlined in section 12.11 
by construction Gj(k) is a better approximation to the exact result in the region of 
reciprocal space close to each of the points {K„}; (in fact other Gj(k) with j^i have 
discontinuities at {K„}j). This suggests that one should construct a new Green's 
function G^(k) which follows each Gj(k) close to each of its corresponding set of 
cluster momenta {K„}i only, smoothly joining up with the other Gj(k) close to 
their set of {K n }i. Such a combined solution G_M-(k) is guaranteed to be causal 
since it will always lie between causal extremes Gj(k). Furthermore, the problem of 
discontinuities is automatically removed since no solution would be followed where it 
has a discontinuity. 

In addition to yielding a new mixed Green's function GM(k) which is fully k- 
dependent with no discontinuities, observe that the method effectively increases the 
number of tiles dividing up the BZ. This leads to a better sampling of the BZ and is 
thus expected to yield meaningful bulk quantities such as density of states. The strict 
criteria that the method must satisfy are outlined in the next section. 

2.3. Formalism: Introduction of continuous k- dependence 

In section \2. 21 it was proposed that a new unique solution (Green's function) may be 
constructed by combining solutions associated with the different choices of {K„}. We 
propose the following criteria that the combined solution Gjw(k) must obey: 

(i) Follow the individual solutions through their associated K„-points (hence yielding 
a meaningful DOS). 

(ii) Be based on a well founded mathematical principle such as the variational 
principle. 
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(iii) Yield a spectral function Agfa) with no discontinuities. 

(iv) Be causal everywhere in k-space. 

(v) Conserve the total number of electron states N. 

(vi) Retain the point-group symmetry of the lattice. 

Here we propose that a general mixed solution G^f (k) satisfying all the above criteria 
may be written in the form 

G M (k) = ^ 7l (k)G l (k) (7) 

i 

where 7i(k) is the mixing parameter of the solution Gj(k) defined by (0J, with 
i = P, AP.. etc. This mixing parameter is generally k-dependent with range 

0<7,(k)<l (8) 

and must satisfy 

E^w = 1 ( 9 ) 

i 

For example, for the case of a 2D square lattice the new Green's function may be 
written as 

G M (k) - 7 p(k)Gp(k) + 7AP (k)G AP (k) 

+ 7p v ap x {k)G PyAPx (k) + y PmA p v (k)G P=cAPy (k) (10) 

where each 7^(k) is to be determined. The new Green's function G^/(k) will obviously 
yield a new set of bulk properties and a unique DOS with its own Fermi energy Ep. 
The spectral function can then be defined in terms of Gm (k) by 

A^(k) = --ImG M (k) (11) 

TT 

Using Gm(1c), a new self-energy Ejvf(k) may also be defined from Dyson's equation 
S M (k)=G - 1 (k)-G M 1 (k) (12) 

In |Appendix A| it is shown that G_w(k) satisfies criteria (iv) and (v) above for any 
set of mixing parameters 7i(k) satisfying JHJl and Next, in section l2~H it is shown 
generally how to choose a specific set of 7i(k) so that the remaining criteria are 
satisfied. 



2-4- Formalism: Choice of 7(k) 

Given the functional requirements described in section ^. 3l that the 7i(k) must satisfy, 
one may narrow down a choice of suitable functions. First, Gjf(k) must follow 
each Gi(k) as close as possible to its associated set of cluster momenta {K n } i . As 
explained in section 12.21 and is clear from figure ^ this enables one to remove all 
discontinuities since no solution will be followed at its corresponding tile boundary 
where the discontinuities occur. Specifically, Gj(k) must equal zero at the tile 
boundaries associated with its {K„}; and so each 7i(k) must be a function that is 
equal to 1 in the tile centres {K„}i and smoothly change to a value of 7i(k) = on the 
tile boundaries. Thus the discontinuities in Gj(k) are removed by cancellation with 
7i (k) . d (k) must also be a periodic function in k-space with the same periodicity as 
{K n }j. Therefore 7i(k) should be expandable in a Fourier series made up of terms 
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from J3J). We also want to base our choice on a mathematical principle that gives a 
unique solution. 

Generally for this system we can define a potential fj, that is proportional to the 
density of states (DOS) 



]TG 4 (k) 74 (k)dk (13) 



For example if a power of the system energy E n is used as the proportionality factor 
one can immediately see that the moments (e.g. band energy) will be directly affected 
by the choice of 7i(k). Hence, the choice of 7»(k) cannot be considered arbitrary, 
one must aim to minimize /x through 7i(k). At any point in k-space there will be an 
associated potential 

^(k)oc- 7i (k). (14) 

between Gjif(k) and Gj(k). Each /i 7i (k) has well defined limits, since at its K„-points 
it should be at a minimum, where the NLCPA coarse-grained Gi(K n ) are considered 
best. Also it will be a maximum on the tile boundary where Gj(K n ) is undefined 
and is at its worst. From section I2"3l 7i(k) is independent of the value G;(k), and 
hence our mixing potential (|14f> will have the same proportionality factor for each 
7i(k) which we can denote by the value r. The total potential of 7»(k) is hence given 

by ' 

/x 7 (k)=r^-7i(k) (15) 

i 

Notice that /i 7 (k) at the K„ points will be solely attributed to the associated Gj(k) 
at that point and vice versa on its tile boundary and we have assumed no functional 
form for 7i(k). With the boundary conditions on how each 7i(k) must behave well 
defined we seek to minimize Ijl5|l . To solve this problem one must notice that (I15|l 
can be mapped directly to the potential of a static membrane |18j described by the 
general solution of the wave equation u(k), related to 7i(k) by 

l/2Vu(k).Vw*(k) oc ^ 7l (k) (16) 

i 

since l/2Vu(k).Vu*(k) to a first order is proportional to that of a potential stored 
in a static membrane ^S]. Hence the sum of 7i(k) is equal to the differential surface 
area of the membrane mapped out by u(k). We wish to minimize the total potential 
^^/^(k) while treating all 7i(k) equally. Minimization of (|16|l gives the well known 
Eulcr (Laplace) equation describing wave motion 

V 2 u(k) = (17) 

that is solved generally for a membrane of the same dimensions as the tile for the 
problem i.e. we only allow waves that can propagate according the boundary conditions 
of 7 i(k). 

Example: 2D square lattice As an example, for the 2D square lattice u has the general 
form 

u(k) = Aexp(i(L x /2)k x ) exp(i(L y /2)k y ) (18) 
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where the origin is defined at the T-point in the first Brillouin Zone, L x and L y are 
the lengths of the cluster (principle vectors) along the x and y directions respectively, 
and A is a complex constant. Substituting l|18|) into l|ltj|) gives 

X>00 oc AA*((L x /2f + (L y /2) 2 ) x 

i 

[cos 2 ((L x /2)k x ) + sin 2 ((L x /2)k x )} x 

[cos 2 ((V2)fc,)+sin 2 ((V2)fc,)] 
= AA*((L x /2) 2 + (L y /2) 2 ) x 

[cos 2 ((L x /2)k x ) cos 2 ({L y /2)k y ) 

+ sin 2 ((L x /2)k ) sin ((Ly/2)k y ) 

+ cos 2 ({L x /2)k x ) sm 2 ({L y /2)ky) 

+ cos 2 '{(Ly I '2)k y ) sin 2 {{L x /2)k x )\ 

= AA*((L x /2) 2 + (L y /2) 2 ) (19) 

where 7i(k) can now be determined by comparison of H19fl and © to yield 

7P (k) = (cos((L x /2)k x )cos{{L y /2)ky)) 2 

7^ P (k) = (sm((L x /2)k x )sm((L y /2)ky)) 2 

j PyAPx {k) = (sm((L x /2)k x )cos((L y /2)ky)) 2 

1p x ap v (k) = (cos((L x /2)k x ) sm({Ly/2)ky)f. (20) 

Each j{ (k) can be seen to uphold their respective boundary conditions given in section 
12.31 Each solution to G^(k) is thus treated equally while removing all discontinues 
and minimizing the potential /i 7 (k), and that of /i given the restrictions on 7i(k). 

Note that each 7^(k) smoothly varies between a value of 1 at {K„}i and at the 
associated tile boundary and has the same periodicity as the tiling. Also, because the 
sum over ji (k) in © is always satisfied even at the tile boundaries one can see that all 
discontinuities are removed because the solution remains continuous and bounded (lies 
in the range between the smallest and largest solutions of Gi (k) given by ©). At the 
k-points where there are no {K„},; the method will take a mix of Gi(k) with weights 
corresponding to the locality of their {K„}i points. It is also worth mentioning the 
limiting cases of N c = 1 and iV c — > oo. When N c — 1 all Gi(k) will equal the same 
CPA value. This is because there will only be one tile (which encompasses the whole 
Brillouin Zone) where 7i(k) always satisfies ©, thus GA/(k) will also equal the CPA 
value. When N c approaches infinity all Gi(k) will approach the exact solution for 
G(k) given by © hence so must Gm (k) because © bounds them within the range 
of Gi (k) . Also note that (|18|) can be extended to 3-dimensions with the inclusion of a 
k 2 term of the same form of the other dimensions, where all mathematical properties 
of 7i(k) shown for the 2D case will hold. 

3. Results 

3.1. Validity of the new method 

The ideal test case for the validity of the new method is the simple ID tight- 
binding model investigated in jSJ. This is because in ID the exact result can be 
obtained numerically. Also, fluctuations are much more significant in ID and so 
detailed structure is expected in the DOS which can be accurately compared with 
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Figure 2. Configurationally-avcraged density of states (DOS) per site as a 
function of energy E for various cluster theories (top to bottom) with cluster 
sizes N c = 2, 4, 6 (W = 1, Vj = ±2) (left to right). 
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both the exact result and other cluster theories such as the embedded cluster method 
(ECM) ^15] and molecular coherent-potential approximation (MCPA) JU] . Indeed, in 
the absence of chemical short-range order the NLCPA results should be very similar 
to those of the ECM and MCPA since the effects of self-consistency are less significant 
in the random case. Note that in ID there are only two NLCPA boundary conditions 
that are possible (periodic and anti-periodic) and the new mixed Green's function is 
given by the simple expression 

G A/ (k) = G P (k) cos 2 ({L x /2)k x ) + G AP (k) sin 2 {{L x /2)k x ) 

= G P (k)cos 2 ((A c a/2)k) + G AP (k) sin 2 ((A c a/2)k) (21) 

Figure 3 in [Sj shows the exact DOS plot for a random A^qB^q alloy of constituent 
site energies Va = +2.0 and Vg = —2.0 (or written more concisely Vi = ±2) with 
hopping parameter W = 1.0, together with the conventional CPA result. Figure[2]in 
this manuscript shows results for the same model obtained using the ECM, MCPA 
(averaged over all cluster sites) and the NLCPA for cluster sizes N c = 2, 4, 6. The 
NLCPA anti-periodic N c = 2 result is the same as the CPA here due to symmetry. 
First note the large difference between the periodic and anti-periodic NLCPA results 
for these small cluster sizes. Significantly, both the periodic and anti-periodic results 
only resemble the ECM and MCPA over particular and different energy regions. 
Indeed, troughs are present at certain energies which are not seen in either the ECM 
or MCPA results or the exact result given in figure 3 of |B| ■ These unphysical features 
are caused by the discontinuities in the Green's function. On the other hand, the 
new mixed NLCPA result reassuringly looks remarkably similar to both the ECM and 
MCPA results for all cluster sizes, with unphysical features removed and troughs and 
peaks in the DOS which can be associated with specific cluster disorder configurations. 
Note that all NLCPA results become increasingly similar as the cluster size increases 
and eventually converge to the same result at the critical cluster size of N c = 12 6 . 

3.2. Fermi surface of a two-dimensional model alloy 

In the 2D case of the NLCPA (and DCA) for systems with Green's functions of the form 
spectral features are unphysical for small cluster sizes due to the coarse graining 
procedure. By having only a small number of tiles any large discrepancy in the self- 
energy will result in a large discontinuity at the tile boundaries and hence a Bloch 
spectral function that remains essentially coarse grained in each tile. Furthermore, 
this behaviour will be most evident near energies that represent constituent ordered 
and clustered alloy states. Bulk quantities may as a result display features that 
overcompensate for short range correlations. For example in the DCA the stable 
d-wave state observed for N c = 4 in the doped repulsive 2D Hubbard model 20 was 
later argued to be unfounded in comparison to exact results [21] . The reason for this 
is that the N c = 4 periodic case exaggerates cluster induced d-wave behaviour by 
the coarse graining procedure. Hence the DCA using periodic boundary conditions 
can only be deemed to correctly predict such features for large cluster sizes [22 ■ It 
is very important from an alloy point of view to improve upon the results available 
from small periodic clusters, since a first principles treatment of large clusters would 
require an unrealistic amount of computational power (at the present time N c < 4 
cluster sizes have been investigated [HI El)- In this section we again backup the validity 
of our method with a DOS comparison, then investigate some of the new interesting 
Fermi surface physics possible with the new implementation of the NLCPA. Also note 
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Figure 3. (colour online) (top) Split-band N c = 4 (W = 0.25, V, = ±1) 
density of states (DOS) using periodic (P), antiperiodic (AP), anti-periodic along 
x-axis (PyAPx), anti-periodic along y-axis(P x AP y ), and mixed (M) solutions, 
calculated from their respective Green's functions. Where the DOS for (P y AP x ) 
and (P x AP y ) are the same by symmetry. To demonstrate the bounding of the 
mixed Green's function in k-space, cross-sections of As(k) at E = — 1 are also 
plotted along the high symmetry axis r — X (centre) and T — M (bottom). 
The mixed solution Ag(k) can be seen to remove discontinuities introduced by 
the original NLCPA coarse graining procedure at tile boundaries e.g 7r/2a. Any 
artefacts of the original coarse graining are removed while correctly incorporating 
the main spectral features of the separate Gi(k). 
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Figure 4. (colour online) A B (k) at E = -1.02 and iV c = 4 (W = 0.25, V* = ±1), 
plotted for the CPA (top-left) and NLCPA-Mixed (top-right) in k-space. Also 
shown is the plot zoomed in to the first quadrant of the Brillouin zone for 
CPA (bottom-left) and NLCPA-mixed (bottom-right) and the anisotropic 
difference A ( ^ PA ('k) — A^(k) (bottom-centre). The colour scale for A B (k) 
represents white as peaks going through yellow and red to troughs as black. The 
colour scale for the anisotropic difference A^ PA (\i) — AM (k) represents positive 
peaks as yellow going through red and blue to negative peaks as black. 



that we have used the concise convention of labelling the Va and Vb atomic potential 
parameters by Vi ± Va since we only use atomic potentials symmetric about E = 
i.e. V B = -V A . 

For this simple alloy model we investigate the first available cluster after N c = 1 
(CPA) i.e. N c — 4 with a nearest neighbour hopping interaction W — 0.25 and cluster 
principle vectors L y = L x — 2a. We maintain the same disorder and constituent 
alloy concentrations as in the previous section for the Id model i.e. a random A50B50 
alloy. The 4-site cluster is the ideal test case since the difference between d (k) should 
be large and so any improvement given by Gm(^-) over Gp(k) should be clear. The 
interaction W has been given the same value as previous work in order to compare 
with the DOS calculated there for the periodic iV c = 32 case. For this lattice all 
solutions for Gi(k) can be considered almost indistinguishable from each other at 
N c = 32 0E], and hence provides a reliable guide to what the exact DOS should 
be. We first examine the DOS (figure |3J for the case of Vi = ±1 where the band has 
split and the nonlocal effects should be large, where we will denote the periodic DOS 
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Figure 5. (colour online) (top) N c = 4 (W = 0.25) density of states (DOS) 
plotted for NLPCA-Mixed, and CPA for varying atomic potentials Vi. Also 
shown (bottom) are the Fermi energy plots of Ag(k) in the first quadrant of the 
Brillouin zone for the half-filled band (E = 0) for the CPA and NLCPA-Mixed, 
where the colour scale represents white as peaks going through yellow and red to 
troughs as black. As Vi is increased a significant change is observed between the 
CPA and NLCPA Fermi surface topology. While the CPA retains a Fermi surface 
reminiscent of a van Hove singularity, the mixed NLCPA produces distinct necking 
features around the X-points. Looking at the DOS we see a small fluctuation in 
the NLCPA-Mixed profile during the transition where the DOS becomes lower 
than the CPA at the Fermi level. 
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by DOSp etc. As one might expect we see that DOSp shows large fluctuations around 
DOSc pa (the same as DO Sap by symmetry) that are unrealistic when compared to 
previous work [2|. Particulary the extra peak formation is much stronger, and results 
in the DOS being far too low at the troughs between peaks with 2 extra semi-bands 
being formed near the main band gap. Indeed in the literature, for the N c = 32 case 
there is never such a large drop below the DOSc pa- Looking now at the new mixed 
result DOSm one sees that the large oscillations associated with DOSp are no longer 
present but the structure is much more like that of the N c = 32 case. Of note is 
how the large troughs of DOSp are filled in: the DOS can be seen to fill in the band 
gap, as well as the central peak becoming prominent over the adjacent structure. In 
summary we see there is a shifted and heightened peak at E « ±1.02 and a second 
peak emerging at E ss ±0.79, that overall corresponds to a clear improvement over 
all N c — 4 solutions separately as well as the CPA. This is more strong evidence that 
the mixing of NLCPA solutions is physical for this model when carried out using the 
method described in this paper. 

The structure of DOSm is a consequence of how Ag(k) behaves in k-space as 
shown in figure|Hl By following each solution according to the weighting function 7i(k) 
the new solution GmO*) removes discontinuities at all tile boundaries (irrespective of 
cluster size) giving a smooth representation of the Bloch spectral function. In figureOH 
cross-sections of the spectral function A^ (k) for an example energy show how the new 
solution is bound and how it is an improvement. Where discontinuities are present, 
the solution switches and uses a supposition of the ones that remain analytic. For 
example Ag (k) at the T-point must be equal to the periodic solution Ag (k) since all 
other solutions have a tile boundary at that point, and so their associated {7i(k)} = 0. 

A comparison can be made for the new spectral function Ag (k) taking the Fermi 
energy at Ep = —1.02 (corresponding to the main DOS peak of the mixed solution) 
with the CPA spectral function for the same energy value, as shown in figure 0] 
Again, some quite remarkable features are apparent, most notably the topological 
change induced by a kink halfway along the X — X direction (not possible in the 
CPA) giving rise to a transition from a circular-like surface to a diamond-like surface 
centred on the T-point, and vice versa for the surface centred on the M-point. This 
is quite dramatic since a Fermi surface (defined spectral function peaks) taken at this 
energy would now have a whole new set of nesting vectors such as the vector along 
X — T — X. Moreover, these 'rounding' and 'squaring' features of Fermi surfaces 
not seen in the CPA calculations have been observed experimentally. For example, 
the work of Wilkinson et al 11 gives a clear example for Gu\^ x Pd x where the 
diamond-like Fermi surface on the (100) plane predicted by the KKR-CPA is actually 
found experimentally to be rounded. Furthermore, smearing and the associated width 
of peaks regains a k-dependence that produces so-called 'hotspots' along the T — X 
directions as shown by dark patches in the difference between the CPA and mixed 
NLCPA. This effectively produces necking in the surface at the X-point in the direction 
X-M. 

In the final set of results shown in figure [5] we examine the Fermi surface for the 
half-filled band with Ep = 0. This is the same for all N c = 4 NLCPA solutions and 
the CPA by symmetry. A Fermi surface is only visible at this energy when there is 
finite DOS and so we examine the properties where there is a breakdown of Virtual 
Crystal (VCA) behaviour to a two-band system. This is achieved by increasing the 
value of the atomic potentials Vi from zero. In the CPA one sees a transition from a 
band with a central peak in the DOS, associated with a van Hove singularity at Ep, 
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to two peaks close to Vi that represent the inclusion of states from both constituent 
metals However, although this was a great breakthrough at the time of the 

discovery of the CPA, no such behaviour is observed in k-space. Sitting in the overlap 
of metal bands one would expect the Fermi level to cross electron states from both 
metals and hence produce a more complicated Fermi surface. However, because the 
CPA has no k-dependence in £ the correction to the VCA only emerges through an 
increased smearing cx Im{Yj(EF = 0)}. There is no topological change in this 
system since i?e{E(£ , F = 0)} = 0. This is because the alloy potentials are equidistant 
from zero and so neither constituent has a preferential weighting. The CPA Fermi 
surface is the same as the VCA Fermi surface at this energy and alloy concentration 
i.e. a straight line joining X-points. With the inclusion of k-dependence in the NLCPA 
a transition to a more complicated form should be possible. However, associated with 
such a transition (near a van Hove singularity) are very large fluctuations in the 
correction to k-space between coarse grained tiles that are not identical. Hence, the 
standard (periodic) solution to the NLCPA does not produce a physical result in this 
regime. However, dramatically, we can now start to investigate such phenomena using 
the new mixed solution to the NLCPA. Results are plotted in figure El for increasing 
Vi and show how a topological transition occurs at Vi w ±0.39 from a CPA-like 
Fermi surface. The Fermi surface splits by introducing new circular surfaces around 
M and T-points and also new necking features around the X-points, in contrast to 
the CPA which only displays increased smearing. This is very important because of 
the new nesting vectors introduced The emergence of this necking around 

X-points is very similar to that observed in recent work for the shape-memory alloy 
Ni^Ales [25, ■ Comparison of first-principles KKR-CPA calculations and experiment 
show the KKR-CPA does not predict a necking feature in the Fermi surface at the 
X-points of the Brillouin zone. The KKR-CPA retains a flat surface associated with 
a van Hove singularity; the topological discrepancy is argued to be the result of the 
Fermi level being close to a van Hove singularity in the DOS. As Vi is increased further 
the surface around M becomes prominent and may be associated with states on the 
band edges of the pure A and B metal states emerging i.e clustering of A's or S's, 
and it is believed X-point features are associated with ordered alloy states. These 
points will be further investigated in future work based on this paper |24| which will 
also incorporate chemical short-range order (which can be introduced by weighting 
the cluster probabilities) and concentration dependence A X B\- X . 

4. Conclusions 

This paper has presented a method for calculating the spectral function within the 
NLCPA. The method is general and can be applied to any (allowed) cluster for any 
lattice in any dimension, and the obtained spectral function is fully k-dependent, 
continuous, and causal. Furthermore, we believe the method to be the correct 
implementation of the NLCPA for cluster sizes below the critical cluster size (see 
introduction) since it yields meaningful results for bulk properties such as the DOS 
and converges systematically as the cluster size increases. 

As an illustrative example, we have applied the method to a 2D tight-binding 
model alloy. Dramatically, a topological transition is clearly observed in the Fermi 
surface as the difference in core potentials is varied in the split-band regime. 
Furthermore, Fermi surface splitting has been observed producing new nesting features 
previously not observed in the conventional CPA. This occurs even in the absence 
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of chemical short-range order (SRO) which could be straightforwardly included by 
appropriately weighting the ensemble of cluster probabilities. 

Clearly, it would be beneficial to apply the method within first-principles 
calculations. Indeed, a derivation and implementation of the spectral function has 
recently been given within the first-principles KKR-NLCPA multiple-scattering theory 
in [25J - However, only the periodic set of cluster momenta {K„} are used and no 
attempt is made to remove the discontinuities. It is our intention to combine the 
method in this paper with the formalism of [21], and investigate the effects of SRO on 
Fermi surface features of systems such as CuPd [TT] , 

Finally, note that the one electron spectral function of correlated electron 
systems also have finite widths. Moreover, they can be measured in angle-resolved 
photoemission experiments (ARPES). As they are due to dynamical fluctuations 
related to the many-electron system they can be described by the DCA. Indeed, the 
DCA has recently been used to interpret such experiments . The method proposed 
here for calculating the spectral function in the NLCPA can be directly transferred 
to the analogous calculation in the DCA and thus it may be useful in this context as 
well as in the study of random alloys. 
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Appendix A. Proof: causality and conservation of states 

It is straightforward to see that the new mixed Green's function Gjif(k) is always 
causal. Provided the set of mixing parameters satisfy JSJ and ©, it follows 
immediately from that Gj\/(k) will always be bound between the individual 
solutions Gi(k) in k-space. Since these individual solutions have an imaginary part 
that is positive definite and hence are causal [2], this implies that Gm(1c) must also 
have an imaginary part which is positive definite and hence is causal. 

In order to prove that the mixing scheme preserves the total number of electron 
states N, we introduce 7i(k) as as a coarse grained quantity in k-space 7i(k) = 7i(K 7 ), 
where {K 7 } are an arbitrary set of fc-points within the BZ that are each centred on a 
coarse grained tile. If we also define M c as the total number of tiles (independent of 
N c ), we have a step function < 7j(K 7 ) < 1 that has M c possible values. For M c = 1 
we have the simplest form of averaging with 7i(k) having no k-dependence, equivalent 
to a weighted average over the DOS. For M c > 1 coarse grained k-dependence is 
introduced to the mixing of G^(k), and M c — > oo representing the introduction of 
complete k-dependent mixing. We wish to show that conservation of electron states 
applies to all values of M c . To do this we integrate Q) over the coarse grained 7i(k) 
for an arbitrary M c to obtain 



where the small k-dependence has been integrated out, leaving a new coarse grained 
version of Gjvr(k), Ga/(K 7 ). Integrating l|A.l|) over all energies — oo < E < oo, we 
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obtain 

N M (K n ) = 5>*(K 7 ) Ni (K n ) , (A.2) 

i 

where Am(K„) and iVj(K n ) are the total number of electron states for the mixed and 
i-th NLCPA solutions respectively for the tile centred at K 7 . Note that this requires 
7(K 7 ) to be energy-independent to keep Nm in terms of N% i.e. we do not want to 
infer any new self-consistency condition on Nm to conserve N apart from the NLCPA. 
Since each solution satisfies the conservation laws separately and electron states are 
evenly distributed across the Brillouin zone, we have JVj(K n ) = N/M c for each i. Now 
(|A.2|) becomes 

i 

Provided © and © are satisfied so that GmO*-) is bounded between the individual 
solutions at every point in k-space, we have for any M c , Nm = N. In particular we 
are free to use the M c — > oo limit to introduce fully k-dependent mixing of G^(k). 
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